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Abstract 

A thermoacoustic engine (TE), which converts thermal energy into acoustic power by the thermoacoustic effect, shows several 
advantages due to the absence of moving parts, such as high reliability and long lifetime associated with reduced manufacturing costs. 
Power output and efficiency are important criteria of the performance of a TE. In order to increase the acoustic power output and ther¬ 
mal efficiency of a Stirling TE, the acoustic power distribution in the engine is studied with the variable load method. It is found that the 
thermal efficiency is independent of the output locations along the engine under the same acoustic power output. Furthermore, when the 
pressure ratio is kept constant at one location along the TE, it is beneficial to increasing the thermal efficiency by exporting more acoustic 
power. With nitrogen of 2.5 MPa as working gas and the pressure ratio at the compliance of 1.20 in the experiments, the acoustic power 
is measured at the compliance and the resonator simultaneously. The maximum power output, thermal efficiency and exergy efficiency 
reach 390.0 W, 11.2% and 16.0%, which are increased by 51.4%, 24.4% and 19.4%, respectively, compared to those with a single R-C 
load with 750 ml reservoir at the compliance. This research will be instructive for increasing the efficiency and making full use of the 
acoustic energy of a TE. 

© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 

Wide attention has been focused on the TE for its poten¬ 
tially high efficiency, simple configuration, reliability and 
environmentally friendly working gas. A TE is an ideal dri¬ 
ver for pulse tube cryocoolers (PTCs) and other thermo¬ 
acoustic refrigerators (TRs). Power output and efficiency 
are two important parameters of a TE. Up to now, most 
attention has been focused only on the acoustic power out¬ 
put [1,2], and the acoustic power output is measured with 
one load each time in most literatures, which is called the 
single output method in this paper. In fact, the interior 
acoustic field of a TE is the basis for power output. There 
is a close relationship between the acoustic power distribu¬ 
tion inside the engine and the acoustic power output, which 
should be paid more attention. 
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A Stirling TE possesses a torus and a branch resonator, 
which is more complicated than the standing wave TE [3,4]. 
It can obtain higher efficiency for its reversible thermo¬ 
acoustic conversion processes. The engine converts thermal 
energy into acoustic power in the thermoacoustic core, 
which comprises an ambient heat exchanger, a regenerator 
and a hot heat exchanger. Then, the acoustic power is 
transferred to the tee where it is divided into two parts. 
One part goes to the resonance tubes, and the other part 
returns along the feedback tubes. Theoretically, acoustic 
power can be exported simultaneously from different loca¬ 
tions along either the resonance tubes or the feedback 
tubes. In this paper, the effects of R-C loads at different 
locations on the acoustic power output and acoustic power 
distribution inside the engine are studied. Then, we propose 
to set two R-C loads along the engine simultaneously to 
increase the acoustic power output and the thermal effi¬ 
ciency, which is called the double output method, to be dis¬ 
tinguished from the single output method [5]. As is verified 
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by the computational and experimental results, the power 
output and the efficiency of the engine with double outputs 
are significantly increased. 

2. Computational results 

2.1. Introduction of the engine and computation method 


Reservoir 



Fig. 1 shows the schematic diagram of the Stirling TE. 
Detailed configuration parameters and performance have 
been introduced in Refs. [4,6]. Acoustic power is exported 
from two locations. Location 1 is at the compliance in the 
torus, and Location 2 is at the start point of the resonator. 
The R-C variable load method is adopted for acoustic 
power computation and measurement [2,7]. Fig. 2 shows 
the schematic diagram of the R-C load, which consists of 
a reservoir, a needle valve and two pressure sensors. In 
the computations and experiments, the volumes of the 
two reservoirs are 500 ml and 750 ml, respectively. 

Since the dimensions of the R-C load are much smaller 
than the acoustic waveguide of the TE, it can be modeled 
as two lumped impedances, an acoustic resistance and an 
acoustic compliance. The compliance impedance can be 
calculated from the reservoir volume V and working condi¬ 
tion. If the dynamic pressure p T is known, the volumetric 
flow U through the needle valve can be determined as 


icoVpr 

yPi n 



where i is the imaginary unit, y is the heat capacity ratio, co 
is the angular frequency and p m is the mean pressure. The 
acoustic power dissipated through the needle valve is de¬ 
fined as 


Fioad — 


oiV\p T \\p E \ sin 6 
2 7Pm 


( 2 ) 


where \p E \ and \p T \ are the pressure amplitudes before and 
after the valve and 6 is the phase difference between p E and 
p T . The resistance of the valve can be calculated by 


Fig. 2. Schematic diagram of acoustic power measurement. p E and p T 
show the locations of the two pressure sensors. 
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Apparently, the acoustic power dissipated through the 
valve is the acoustic power output of the TE. 

The computation is performed by using DeltaE (design 
environment for low amplitude thermoacoustic engines), 
which was developed by W.C. Ward and G.W. Swift. 
Based on linear thermoacoustics theory, DeltaE includes 
modules for heat exchangers, regenerators, ducts, pistons 
etc., so that a complete Stirling TE can be modeled [8,9]. 
We model the self made TE (shown in Fig. 1) as 64 mod¬ 
ules connecting in turn, including the torus and resonance 
tubes. Viscous dissipation and nonlinear dissipative effects, 
such as dissipation associated with bends and tapered 
ducts, are considered. Nitrogen of 2.5 MPa is used as the 
working gas. 

2.2. Effect of output location with single output method 

The engine shown in Fig. 1 is a 1/4 wavelength Stirling 
TE. Its pressure ratio becomes maximal just above the 
main cooler, decreases gradually anticlockwise along the 
torus and reaches the pressure node before the tapered res¬ 
onance tube [10]. So, the acoustic power output should be 
different at different locations due to the acoustic field dis¬ 
tribution. In the computation, the pressure ratio at the 
compliance is kept at 1.20. 




9 Resonator; 


Fig. 1. Schematic diagram of the thermoacoustic engine. 
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Fig. 3 shows the acoustic power output dependence on 
resistance with the single output method at different loca¬ 
tions. All the computational output curves have similar 
trends. The curves peak when the resistance is in the range 
of 3.0-5.0 x 10 7 N s/m 5 . The compliance impedances of the 
500 ml and 750 ml reservoirs are 4.55 x 10 7 N s/m 5 and 
3.03 x 10 7 N s/m 5 , respectively, which indicates that the 
acoustic power output reaches its maximum when the resis¬ 
tance and the compliance impedance of the R-C load are 
equal. The computational results also show that the engine 
can export more acoustic power with a reservoir of larger 
volume and a valve of larger flow coefficient (i.e. smaller 
resistance). The acoustic power outputs at the compliance 
are much larger than those at the start point of the reso¬ 
nance tubes, which is due to the acoustic field distribution. 
When the resistance of the R-C load is larger than 
1.0 x 10 8 N s/m 5 , it has a dominant influence on the output 
characteristics compared to that of the compliance imped¬ 
ance. So, the output curves with different reservoirs at the 
same location are almost superposed. 

Fig. 4 shows the acoustic power distribution in the TE 
with the R-C load with 500 ml reservoir at the start point 
of the resonance tubes. When the resistance is less than 
1.0 x 10 5 N s/m 5 or larger than 1.0 x 10 10 N s/m 5 , about 
330 W acoustic power is transferred to the resonator and 
dissipated to maintain the interior acoustic field. About 
500 W acoustic power is generated by the thermoacoustic 
core, while the acoustic power returned to the ambient 
end of the regenerator is about 875 W, which indicates that 
more acoustic power can be exported from the torus. Fig. 5 
shows the acoustic power distribution with the R-C load 
with 750 ml reservoir at the compliance in the torus. More 
acoustic power is returned at the tee, and more acoustic 
power is exported as expected. The more important result 
is that the TE generates much more power in the thermoa¬ 
coustic core when the power output curve rises. As is 
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Fig. 4. Acoustic power distribution in the engine with the R-C load with 
500 ml reservoir at the start point of the resonance tubes. 
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Fig. 5. Acoustic power distribution in the engine with the R-C load with 
750 ml reservoir at the compliance in the torus. 



Fig. 3. Comparison of single power outputs at different locations. Solid, 
dash, short dash and dot curves are computational results. Scattered 
points are experimental results. 500 ml and 750 ml denote the related 
reservoirs of 500 ml and 750 ml for conciseness. 


shown, when the power output peaks, the feedback acous¬ 
tic power at the main cooler reaches its minimum, while the 
power generated in the core peaks. Comparison with Fig. 4 
indicates that the regenerator works efficiently here. It is 
also shown that the acoustic power transferred to the reso¬ 
nator is kept nearly constant when the variable R-C load is 
set at the compliance. 

In order to make clear where acoustic power can be 
exported more efficiently, the dependences of the thermal 
efficiency and exergy efficiency on acoustic power output 
at different output locations are shown in Fig. 6. The defi¬ 
nitions of thermal efficiency and exergy efficiency of a TE 
are 



( 4 ) 

( 5 ) 


























D. Sun et al. / Energy Conversion and Management 49 (2008) 1265-1270 


1268 


o3 


<D 

4=! 

H 


30 


25 - 


o 
C 

<D 

o 

4h 
<D 

6D 

S-H 
QJ 

X 

<L> 

a 10 

o3 


20 - 


15 - 


5 - 


0 


Computational thermal efficiency 

.750 ml at location 1 

- 500 ml at location 2 

Computational exergy efficiency 

- 750 ml at location 1 

- 500 ml at location 2 


A 


A* 



A • 


Experimental thermal efficiency 
^ 750 ml at location 1 

• 500 ml at location 2 

Experimental exergy efficiency 
A 750 ml at location 1 
O 500 ml at location 2 

_i_i_i_i_i_i_i_i_i_ 


0 50 100 150 200 250 300 350 400 450 500 

Acoustic power (W) 

Fig. 6. Thermal and exergy efficiency dependence on acoustic power 
output. Solid, dash, dot and dash dot curves are computational results. 
Scattered points are experimental results. 500 ml and 750 ml denote the 
related reservoirs of 500 ml and 750 ml for conciseness. 


where W is the net acoustic power output to the R-C load, 
Q is the heating power, r 0 is ambient temperature and T n 
is the wall temperature of the heater. The acoustic energy 
stored in the interior acoustic field is dissipated by the 
TE itself. When calculating the thermal and exergy efficien¬ 
cies, we only take the acoustic power exported by the R-C 
loads into account, which is different from Ref. [11]. If 
there is no net acoustic power output, the efficiency of 
the TE is zero. Two conclusions are drawn from the com¬ 
parison of the curves. One is that the thermal efficiency and 
exergy efficiency are independent of power output location 
under the same acoustic power output. Although the 
acoustic power outputs at different locations are different, 
even with the same R-C load (shown in Fig. 3), the ther¬ 
moacoustic conversion efficiencies are at the same level so 
long as the acoustic power outputs are equal. The other 
is that the thermal efficiency and exergy efficiency increase 
with acoustic power output. For the pressure ratio at the 
compliance of the TE remaining at 1.20, the power dissi¬ 
pated in the whole engine changes little during the whole 
operation, which means that more acoustic power output 
leads to lower average dissipation. Comparing Fig. 6 with 
Fig. 3, we can find that much more acoustic power output 
and higher efficiency can be obtained by exporting acoustic 
power at the compliance. 

2.3. Performance with double output method 

Inspired by the above results, we then use the so called 
double output method to export more acoustic power with 
two R-C loads simultaneously (as shown in Fig. 1). In 
order to obtain the trends of acoustic power output with 
two loads, the valve at the resonator (Focation 2 in 
Fig. 1) is fully opened first (by decreasing the resistance 
in the computation), and the valve at the compliance 


(Focation 1 in Fig. 1) is gradually opened. In the computa¬ 
tion, the pressure ratio at the compliance is also kept at 
1 . 20 . 

Fig. 7 shows the comparison of power outputs between 
the single and double output methods. All the computa¬ 
tional curves peak at particular resistances. On the right 
side of the peak, the total acoustic power output increases 
with the opening of the valve at the compliance. The total 
acoustic power with the double output method is approxi¬ 
mately the sum of the powers with the single output 
method, which implies that the outputs at different loca¬ 
tions have little influence on each other. During the varia¬ 
tion of the resistance at Focation 1, the power output with 
the 500 ml reservoir at Focation 2 varies little. With the 
750 ml reservoir at Focation 1 and the 500 ml reservoir at 
Focation 2, the total acoustic power output reaches its 
maximum. When the resistance and compliance imped¬ 
ances are nearly equal, the maximal acoustic power output 
is about 550 W, while the maximal power output is 430 W 
with the single output method. 

Fig. 8 shows the comparison of the thermal efficiencies 
between the single and double output methods. Keeping 
the resistance the same, the thermal efficiency with the dou¬ 
ble output method is much higher than that with the single 
output method. The highest computational thermal effi¬ 
ciency reaches 23.5% when the resistance is about 
2.9 x 10 7 N s/m 5 , which is increased by 18.7% more than 
the single output thermal efficiency of 19.8%. Comparison 
of the exergy efficiencies is shown in Fig. 9. The maximal 
exergy efficiency with the single output method is 25.7%, 
while it reaches about 30.0% with the double output 
method, increased by 15.6%. 



Fig. 7. Comparison of power outputs between single and double output 
methods. Solid, dash and dot curves are computational results. Scattered 
points are experimental results. 500 ml and 750 ml denote the related 
reservoirs of 500 ml and 750 ml for conciseness. The powers denoted by 
500 ml and 750 ml show the individually exported acoustic powers, 
respectively, with double output method. 
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Fig. 8. Comparison of thermal efficiencies between single and double 
output methods. Solid and dash curves are computational results. 
Scattered points are experimental results. 500 ml and 750 ml denote the 
related reservoirs of 500 ml and 750 ml for conciseness. 



Fig. 9. Comparison of exergy efficiencies between single and double 
output methods. Solid and dash curves are computational results. 
Scattered points are experimental results. 500 ml and 750 ml denote the 
related reservoirs of 500 ml and 750 ml for conciseness. 

3. Experimental verification 

In order to verify the simulation, a number of experi¬ 
ments were performed. The engine and measurement set 
up are shown in Figs. 1 and 2. Four calibrated pressure 
sensors, labeled pi, p 2 , P 3 and p 4 , are arranged in the sys- 

Table 1 

Comparison of experimental results 


tern. Two needle valves (Swagelok®) with maximal flow 
coefficients of 0.73 are applied. In order to keep the gas 
properties uniform, cooling water is used to cool the tubes 
and valves. Nitrogen of 2.5 MPa is used as the working gas. 
The pressure data acquisition system comprises the 
pressure sensors, data acquisition clip, computer and self 
developed program based on Lab view 6.1 by National 
Instruments Inc. (NI) [12]. The pressure sensors are of 
the linear silicon piezoelectric type and were supplied by 
Infineon Technologies, Germany. The data acquisition clip 
is NI PC-1200, 12 bits in precision and 10 k S/s in sampling 
frequency. Calibrated NiCr-NiSi thermocouples, digital 
multimeter (Prema 5017) together with a digital scanner 
are used to measure the temperatures. Heating power is 
adjusted by the charging voltage to the heaters and mea¬ 
sured by a digital dynamometer. 

As in the computation, the pressure ratio at p\ is stabi¬ 
lized at 1.20 by adjusting the heating power. The experimen¬ 
tal data are scattered in Figs. 3, 6-8 and 9, respectively, and 
Table 1 shows the maximal values under different 
conditions. 

It is shown in Fig. 3 that the experimental acoustic 
power output decreases with the increase of resistance 
impedance since the resistances of the needle valves are 
on the right side of the peak. In the experimental range, 
the computational and experimental results are in good 
agreement. As shown in Table 1, with the 500 ml reservoir 
at the compliance, the maximal acoustic power output of 
239.4 W is obtained, while it is only 126.2 W at the resona¬ 
tor. With the 750 ml reservoir at the compliance, the max¬ 
imal acoustic power output of 257.6 W is obtained, while it 
is 125.2 W at the resonator. As shown in Fig. 6, the thermal 
efficiency and exergy efficiency are also increased greatly by 
arranging the R-C loads at the compliance by exporting 
more acoustic power. Resembling the computation curves, 
the scattered experimental points of thermal efficiency and 
exergy efficiency are also superposed in their common 
range, respectively. We did not take serious measures to 
thermally insulate the heat exchangers, the regenerator 
and the thermal buffer tube in the experiments, and this 
heat loss is not included in the computation. Therefore, 
there is a relatively large deviation between the computa¬ 
tional and experimental results. 

As shown in Fig. 7 and Table 1, the maximal total 
acoustic power output with the 750 ml reservoir at Loca¬ 
tion 1 and the 500 ml reservoir at Location 2 is 390.0 W, 
which is increased by 51.4% more than that with the single 



Single output 




Double output 

A* 

500 ml reservoir 


750 ml reservoir 


Location 1 

Location 2 

Location 1 

Location 2 

Maximal power (W) 

239.4 

126.2 

257.6 

125.2 

390.0 

Thermal efficiency 

8.3% 

4.8% 

9.0% 

4.7% 

11.2% 

Exergy efficiency 

12.3% 

7.5% 

13.4% 

7.3% 

16.0% 


A* denotes the condition with 750 ml reservoir at Location 1 and 500 ml reservoir at Location 2. 
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output method. During the adjustment of the valve at loca¬ 
tion 1, the power output with the 500 ml reservoir at Loca¬ 
tion 2 varies little. 

In Fig. 8, the highest thermal efficiency with the double 
output method reaches 11.2% when the resistance is about 
1.0 x 10 8 N s/m 5 , which is increased by 24.4% more than 
the single output thermal efficiency of 9.0%. The relatively 
large deviation between the computational and experimen¬ 
tal results is also caused by the poor thermal insulation. It 
is shown in Fig. 9 and Table 1 that the exergy efficiency is 
increased remarkably too. The maximal exergy efficiency 
with the single output method is 13.4%, while it reaches 
16.0% with the double output method, increased by 19.4%. 

4. Conclusion 

Computations and experiments show that the power out¬ 
put locations have not much influence on the thermal effi¬ 
ciency and exergy efficiency under the same power output. 
Because of the acoustic field distribution, more acoustic 
power can be exported from the locations with higher pres¬ 
sure ratio (or pressure amplitude) with the same R-C load. 
When the pressure ratio at one location is kept constant, the 
acoustic power dissipated inside the engine varies little. So, 
more acoustic power output leads to higher efficiency. 

Compared to mechanical compressors, a TE has a large 
scale interior acoustic field, which makes it more conve¬ 
nient to export acoustic power from more than one loca¬ 
tion along its waveguide simultaneously. For simplicity 
and convenience, we proposed to use two simultaneous 
outputs to show the performance improvements of a Stir¬ 
ling TE in this paper. More acoustic power and higher effi¬ 
ciency have been obtained. From the viewpoint of 
applications, it will be convenient and efficient for a TE 
to drive several TRs at different locations simultaneously. 
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